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The problem of crossflow receptivity is considered in the context of a canonical 3D 
boundary layer (viz., the swept Hiemenz boundary layer) and a swept airfoil used recently 
in the SWIFT flight experiment performed at Texas A&M University. First, Hiemenz flow 
is used to analyze localized receptivity due to a spanwise periodic array of small amplitude 
roughness elements, with the goal of quantifying the effects of array size and location. 
Excitation of crossflow modes via nonlocalized but deterministic distribution of surface non- 
uniformity is also considered and contrasted with roughness induced acoustic excitation of 
Tollmien-Schlichting waves. Finally, roughness measurements on the SWIFT model are 
used to model the effects of random, spatially distributed roughness of sufficiently small 
amplitude with the eventual goal of enabling predictions of initial crossflow disturbance 
amplitudes as functions of surface roughness parameters. 


I. Introduction 

Surface roughness characteristics are known to exert profound influence on the onset of transition in 
swept wing boundary layers. A striking example of surface quality effects on transition is documented in 
Ref. 1, wherein transition is noted to move toward the trailing edge from approximately 25 to 30 percent 
chord to 80 percent chord when the painted leading edge was polished to reduce the roughness amplitude 
from 1.0 fi m root mean square (RMS) (3.8 fi m average peak-to-peak) to 0.3 fi m RMS (2.2 /im average 
peak-to-peak). Had linear stability based N-factor methods been applied to predict the onset of transition 
in this case, the correlating N-factor would have increased from 8 to 16. Similarly large N-factor values have 
been reported in the context of earlier investigations of swept wing boundary layers at ONERA. 1 2 Because 
the small amplitude surface roughness primarily seeds the stationary crossflow instabilities that eventually 
lead to transition, prediction uncertainties of this magnitude can only be minimized by accounting for the 
receptivity phase in the transition prediction process. Inclusion of receptivity is particularly critical in the 
context of laminar flow control via distributed roughness elements (DREs), 3 since the delay in transition is 
achieved via controlled seeding of less unstable crossflow modes and the application of DRE heights outside 
of an optimal range may result in inadequate control at one end (i.e. , too small heights) and premature 
transition (or over-control) at the other end. The latter possibility is particularly likely for high Reynolds 
number configurations as the linear amplification ratios achieve substantial magnitudes across a broad range 
of stationary crossflow wavelengths. 4 

Roughness induced excitation of crossflow instabilities was examined in Ref. 5, wherein preferential ex- 
citation of stationary crossflow vortices was confirmed on the basis of a straightforward extension of the 
Goldstein 6 - Ruban 7 theory. Potentially substantial effects of surface curvature and mean-flow nonparal- 
lelism were pointed out in Ref. 8 and the effects of nonlinear roughness height were recently addressed in Ref. 
9. References 8, 10 19 describe the various applicable prediction methodologies (ranging from the inhomoge- 
neous form of linear stability equations, parabolized stability equations, and linearized or full Navier-Stokes 
equations) and their application to a broader class of flows. Much of this work, however, pertains to a single, 
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spanwise periodic array of roughness elements, while studies reporting the receptivity to spatially nonlocal- 
ized roughness distributions are less common. (See Ref. 20 for a theoretical study of nonlocalized roughness 
effects or Ref. 21 for an experimental investigation on the subject.) For small amplitude roughness, at least, 
the theoretical framework for the crossflow receptivity to distributed surface roughness is similar to that 
of the acoustic excitation of Tollmien-Schlichting waves via nonlocalized roughness. 10 Yet, this framework 
needs to be applied with due consideration to the effect of surface curvature and mean boundary layer growth 
on both receptivity and stability characteristics of stationary crossflow modes. 

An integrated approach combining receptivity with linear and nonlinear growth of the stationary crossflow 
instability was demonstrated in Ref. 19, using parabolized stability equations (PSE) and adjoint PSE as 
implemented in the Langley Stability & Transition Analysis Code (LASTRAC). This combination of tools 
enables efficient and accurate prediction of disturbance evolution in many cases of interest. However, when 
the disturbance wavelength of interest changes rapidly in space, the basic parabolizing assumption becomes 
suspect. There, the full elliptic nature of the flow must be included to achieve an accurate disturbance 
prediction. Harmonic Linear Navier-Stokes (HLNS) tools have been shown to be successful in complementing 
existing PSE tools in these flow environments. (See the work of Streett 13 for a detailed discussion of the 
HLNS approach.) The increased fidelity of the HLNS approach does have its drawbacks: most notably the 
increased computational costs, which can exceed comparable PSE simulation times by at least an order of 
magnitude. Thus, the HLNS approach is relatively less suitable for performing extensive design investigations 
or optimization studies. 

This paper addresses two fundamental challenges. First, a discrete adjoint formulation is introduced to 
complement an existing Harmonic Linear Navier-Stokes capability. Profound improvements in efficiency are 
achieved with the new adjoint approach as shown in numerous examples. Furthermore, it is shown that 
extensive design and optimization studies are now straightforward using the higher fidelity HLNS equations. 
With the new tool in hand, we consider both localized and distributed receptivity mechanisms for crossflow 
vortex instability. Two different swept-airfoil configurations are examined as part of the parameter study. 
First, we consider a canonical 3D boundary layer (viz., the swept Hiemenz boundary layer), followed by 
a realistic swept airfoil geometry employed during a recent SWIFT experiment performed at Texas A&M 
University. 1,22 Localized receptivity due to a spanwise periodic array of small amplitude roughness elements 
is first analyzed using Hiemenz flow, with the goal of quantifying the effects of array location and element 
size. Next, the excitation of crossflow modes via a nonlocalized but deterministic distribution of surface non- 
uniformity is considered, and contrasted with the spanwise periodic roughness elements as well as roughness 
induced acoustic excitation of Tollmien-Schlichting waves. Finally, the SWIFT experimental flow conditions 
are used to study the effects of random, spatially distributed roughness (of sufficiently small amplitude). 
The eventual goal is to predict the initial crossflow amplitudes via a reduced order model of the roughness 
geometry, using only a small number of amplitude and texture parameters to represent the surface roughness. 
It is envisioned that this study will help identify the most relevant parameters and, hence, contribute to the 
potential development of an appropriate database pertaining to the natural roughness of typical aerodynamic 
surfaces. 17 

The paper is organized as follows. Section II gives a general description of the SWIFT flight experi- 
ment. Section III (a) summarizes the existing HLNS formulation originally developed by Streett, 13 while 
Section 111(b) introduces a complementary adjoint formulation. Section IV includes a brief description of 
the theoretical aspects of receptivity. Section V presents results. First, a detailed study of the well known 
incompressible swept Hiemenz flow is presented, including extensive validation of the new adjoint approach. 
Next, the adjoint approach is utilized to compare and contrast the receptivity of the boundary layer to 
spanwise periodic roughness elements, and deterministic nonlocalized roughness distributions. Finally, a 
detailed study of the receptivity over the SWIFT airfoil configuration is presented, beginning first with a 
description of the surface characteristics of the wing. Next, a parametric study of optimal roughness element 
size and location is presented. Finally, the effects of random spatially distributed roughness of sufficiently 
small amplitude is presented. Section VI summarizes the work. 

II. Experimental Setup 

The flow configuration of interest in the present paper corresponds to the SWIFT flight experiment 
conducted at Texas A&M University. 1,22 A 30-degree swept airfoil test article was mounted vertically 
beneath the wing of a Cessna 0-2A Skymaster aircraft and flown at chord Reynolds numbers of 7 to 8 
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million, angle of attack olaoA — —4.7°, and airspeed 180.0 KTAS (knots true air speed), to investigate the 
feasibility of discrete roughness elements (DRE) technology for Sensor Craft applications. The goal was to 
achieve a minimum of 60% laminar flow over the suction side of the model by appropriately placed DREs. 

For the sake of clarity, a single DRE is simply a bump (or, alternatively a dimple) on the surface of 
the airfoil, e.g., in the form of a shallow cylindrical disk (a hockey puck). A spanwise, periodic array of 
such elements is constructed by aligning multiple elements parallel to the spanwise direction, separated by 
a distance d. Thus, the array is periodic with length d. Spanwise-periodic DREs provide a passive means 
of controlling the crossflow instability responsible for transition over a swept wing without the need for an 
active device (e.g. suction). Indeed, Saric et al. 3 demonstrated that roughness elements on the scale of 
micro-meters were effective in moving the transition location. The method for spanwise spacing of DREs 
is outlined in the work of Reibert et ah, 23 whereby appropriately forced modes distort the mean flow to 
inhibit the naturally occurring most unstable crossflow vortices, thus delaying transition. Conversely, DREs 
with different placement and character can be used to excite the unstable modes thus forcing transition and 
achieving a large chordwise extent of turbulent flow. 

Linear stability theory applied to the SWIFT configuration identified the 4.5mm wavelength crossflow 
mode as being extremely unstable, whereby a DRE spanwise spacing corresponding to the more stable 
2.25mm mode could be used for transition control via DREs. The placement of the DREs is typically near 
the neutral point of the control (2.25mm) mode; i.e. between 1% and 2% chord for the inboard and outboard 
stations, respectively. 

For the work described in this paper, the mean flow is computed with a boundary layer solver based on 
the infinite span approximation, using the surface pressure distribution derived from the CFD solution. 22 
Spanwise variations in the computed pressure distribution over the mid-span segment of the suction surface 
were relatively small, so that the region between the rows of static pressure ports at 13 inch and 29 inch 
span locations outboard of the mounting plane are assumed to approximate (locally) infinite span conditions. 
Transition measurements, both without and with the use of DREs, were reported in Ref. 22 High-fidelity 
computational analysis of the DRE based control was reported in Ref. 24. 


III. Analysis Codes 

Many of the receptivity results presented herein are obtained using a harmonic variant of the linearized 
incompressible Navier- Stokes equations. The basic methodology and algorithm was reported in Ref. 13 An 
adjoint formulation has been developed to complement and expedite the original algorithm. A brief outline 
of each approach is now presented. (See Refs. 13 or 25 for further details.) 


III. A. Incompressible Harmonic Linearized Navier-Stokes 

The flows of interest in this study (i.e. M ^ < 0.25) are well represented using the incompressible Navier- 
Stokes (NS) equations. The derivation of the harmonic , linearized incompressible NS equations begins with 
an expansion of the flow around a base state: 


u (x,y,z,t) = u (x,y,z) + e u (x, y, z, t) . (1) 

where U = [P, £/, V, W] T , U = [P, [7, V, W] T , u = [p, u, v,w] T and the components \p,u,v,w\ 
are the incompressible pseudo pressure- and (x,y : z) velocity-perturbations in the x, y and 2 coordinate 
directions, respectively. By construction, the base state U satisfies the steady incompressible NS equations, 
or a suitable approximate form; e.g., the boundary- layer equations or the Falkner-Skan- Cooke similarity 
equations. Substituting equation (1) into the nonlinear incompressible NS equations, and collecting terms 
of order e, yields a linear set of disturbance equations in the variables u (x, p, 2, t; U ). a 

The linear disturbance equation for the variables u (x, p, z, t\ U ) is further simplified by the ansatz 


d u 

dt 


—iuj u 


d u 

dz 


+ i/3 u . 


(2) 


The first ansatz limits solutions to a single harmonic, specified by the frequency ca, while the second assumes 
spanwise periodicity in z at a specified wavenumber /3. Both approximations are well suited for the infinite 


a It is well known that a disturbance equation of this form is suitable for predicting the initial stages of transition, especially 
when amplitudes are orders of magnitude smaller than the background base flow U . 
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span swept wing flows studied herein. The time-harmonic and spanwise periodicity assumptions reduce 
the 4-dimensional original formulation for the variables u (#, y, 2 , t; U ), to one involving only two spatial 
dimensions: for the variables u (x,y;u,/3, U ). 

d(x,y;u>, (3, U) u + a x (x,y;/3, U ) ^- + a y (x,y;/3, U ) ^- = 
b xx (x,y, U) + b xy (x,y; U ) + b yy (x,y ; U) 0- 

Remark: While equation (3) is expressed in terms of the independent variables (x,y), it could equally well 
represent an orthogonal coordinate frame (aq,^). Here, the independent variables (x,y) correspond to the 
orthogonal, streamwise and wall-normal directions. 

Boundary conditions are a sensitive aspect of this formulation contributing both to physical well-posedness 
(accuracy) as well as numerical conditioning. Typically, Dirichlet and/or Neumann conditions are specified on 
the inflow, wall and farfield boundaries. The specification of the down- stream boundary condition requires a 
different approach, however. A buffer-domain technique, 26 in which the momentum equations are smoothly 
parabolized to convect disturbances through the outflow boundary, was found to be robust and reliable 
alternative to conventional Dirichlet or Neumann conditions. 

Solution of the partial differential equations (PDEs) given in equation (3) requires the spatial derivatives 
to be discretized, resulting in a linear matrix equation of the form 

Ahlns w = f 

where w is the solution vector for all discretely orthogonal points (yj,Xi) in the field, and the right-hand 
side vector f contains boundary conditions (e.g. wall blowing/suction, or free-stream disturbances). 

The structure of the matrix Ahlns depends on the form of the discrete spatial derivative operators. In 
the present study, the streamwise direction Xi is discretized using a variable stepsize fourth-order centered 
difference technique, while the wall normal direction yj uses a mapped Chebyshev collocation technique. Grid 
refinement studies were performed in all cases, and established that all results were accurate to approximately 
three significant digits. The maximum tensor product mesh needed was of dimension 4, 000 x 71, while most 
required significantly fewer points (i.e. 1,500 x 50). 

This discretization results in a block pentadiagonal matrix if a lexigraphic (j, i) ordering is used. The 
block size A/ is given by = 4 * ny where ny is the number of yj Chebyshev collocation points. A direct 
solver was used in all cases to invert the linear system. The single-processor CPU times (Intel Xeon operating 
at 3.0 GHz) were in the range 1 — 10 minutes depending on the particular grid. Inversion times scaled linearly 
with the streamwise dimension aq and cubically with the wall-normal yj dimension. The problems considered 
in this work were generally too small to benefit appreciably from multi-threaded parallel approaches (e.g. 
OpenMP BL AS2 and BLAS3 operations) . A modest improvement was achieved by splitting the streamwise 
direction into multiple blocks, coupled using a Schur-Complement approach. This technique enables the use 
of multiple cores on a dual socket quad-core workstation and leads to a reduction in wall-clock time by a 
factor 2 — 3 when using all 8 cores. 

Summarizing, simulations using the quasi-3D harmonic equations are easily preformed on commodity 
hardware ( $1000), thus allowing “routine” engineering calculations in many setting even without any adjoint 
formulation. 

III.B. Adjoint formulation 

Although a continuous adjoint approach was adopted in previous PSE work, 25 herein a discrete adjoint 
formulation is implemented to expedite the solution of equation (3). A brief summary of the contributing 
factors that led to the choice of the discrete approach is now presented. 

III. B. 1 . Background 

There is a long history of the use of adjoint equations in fluid dynamics. Pironeau, 27 29 first used the 
adjoint approach for design studies followed shortly thereafter by design optimization work in CFD (po- 
tential/Euler /NS) by Jameson. 30 32 The breadth and complexity of applications and maturity of the field 
progressed dramatically over the next decade. 33-39 
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Adjoint methods are typically divided into two different categories: discrete and continuous adjoint 
methods, with distinct advantages to each approach. Consider the linear differential equation C subject to 
well-posed boundary data 23, b 

= fi |n ; Bu = f b | 9n (4) 

and an output cost functional J of interest 

J = (g,u)n + (h,Cu)dn • (5) 

Here, the nomenclature (., .) denotes any standard inner product, while and dQ denote action on the 
volume and boundary, respectively. 

The continuous approach seeks a dual formulation 

£> = gi |n ; B'V = h b Ian (6) 

where the adjoint equation £* and adjoint variable ^ satisfy the integral constraint 

(£*Vh u)n + (23>, Cu) dn = (ip, £u)n + (C*ip, Bu) dn (7) 

which is obtained from the original problem using integration by parts (and demanding sufficient smoothness 
of all quantities that are differentiated). With the dual formulation, the equivalent output cost function 
becomes 

J = (Vh fi)n + (C*Vh /&)ao • (8) 

The advantage of the adjoint formulation can be profound when the problem of interest requires multiple 
solutions of equation (4) subject to the datasets [// , f b ; k = 1, • • •]. Rather than one solution with each 
dataset k = 1, • • •, only a single solution of the dual formulation given in equation (6) is required, followed 
by repeated evaluation of the dual cost function given by equation (8). 

Remark. In the continuous approach the adjoint equations and boundary conditions have a precise 
physical interpretation, and are amenable to mathematical analysis. For example, Giles et al. 40 derive 
analytic well-posed adjoint boundary conditions for the Euler equations, and analytically study the behavior 
of these equations in the vicinity of sonic points, a more difficult proposition when using a discrete approach. 

Remark. Note that in the continuous approach the discretization of the equations is explicitly performed 
on the dual equations given in (6) and (8). The discretized continuous adjoint equations (i.e. the matrix) is 
not discretely adjoint to anything. 

In the discrete adjoint approach, the governing equations given in (3) are first discretized, the result being 
a matrix problem of the form AU = F. A mapping is then sought to transfer the original problem over to 
an equivalent dual problem. The discrete adjoint is derived as follows. Given a discrete inner product (., .) 
defined by the vector product 


(£>»?) = 

= fri 


then the original formulation 



AU = F -> 

I = (G, U) 

(9) 

can be expressed using an equivalent dual formulation 



a t * = g -*• 

I = (1>.F). 

(10) 

A proof of the equivalence of the two problems follows 

immediately 


(*,F) = (¥,AU) = 1 

:A t *,U) = (G, U). 

(11) 


Note that equation (10) provides the discrete counterpart of the discretized versions of equations (6) and (8). 
They are not precisely equivalent because the continuous operations of differentiation and integration by parts 
are not in general exactly satisfied at the discrete level. Thus, results obtained by using the two approaches 
may not be identical, the differences being at the level of the truncation error in the discretizations. 

Remark. The physical significance of the discrete equations is more difficult to interpret. At the very 
least however, the adjoint solution provides the sensitivity /influence of an arbitrary source term F on the 
cost functional I of interest. 

b See Giles et. al. 36,37 for further details. 
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III.B.2. Adjoint Boundary Conditions 

The decision to implement a discrete adjoint formulation was strongly motivated by three factors: 

• Boundary conditions. A problematic element of the continuous adjoint approach is resolving the com- 
plexities brought about by the imposition of buffer domain boundary conditions. Despite the avail- 
ability of well-posed Dirichlet and Neumann adjoint boundary conditions for the continuous problem 
(see Giles 40 ), no such conditions exist for the non-reflecting boundary conditions used at the outflow 
plane. Furthermore, the buffer domain approach requires modification of the governing equations in 
the vicinity of the boundary. All these factors must be considered in the derivation of suitable adjoint 
boundary conditions, not an insurmountable task, but subtle nevertheless. 

• Full consistency. The inviscid characteristics terms in the continuous adjoint equations are of opposite 
sign compared with the HLNS equations. Thus, the adjoint outflow boundary coincides with the inflow 
boundary of the continuous equations. Thus, by symmetry the adjoint equations require a nonreflecting 
boundary condition at the upstream boundary, including a modification of the governing equations near 
that boundary. This introduces an inconsistency between the primal and dual equations at each end 
of the domain, unless precise boundary conditions were available. 

• Coding simplicity. The discrete approach only required building the primal matrix Ahlns followed by 
its transposition. Everything is accomplished strictly at the level of linear algebra. Furthermore, the 
primal and dual problems are completely consistent because the same governing equations (including 
buffer domain) are used for both. 

I1I.B.S. Cost Functionals 

The solution of the incompressible Harmonic Linear Navier- Stokes equations in cases dominated by cross- 
flow instability mechanisms is characterized by a growing streamwise disturbance. The amplitude of the 
disturbance (at a fixed location) is a direct measure of the efficacy of various forcing mechanism F. Cost 
functionals that are synonymous with the local nature of the disturbance, provide effective sensitivity mea- 
sures. Indeed, recall that the adjoint solution provides the sensitivity /influence of an arbitrary source term 
F on the cost functional I of interest. Several cost functionals were tested, and all were found to give 
qualitatively similar results. Two were thoroughly tested and then used throughout all studies. 

The first is the integral of the disturbance “energy” at a fixed streamwise location 

rymax 

I u = [u(y,Xi) + v{y,Xi) + w{y,Xi)\dy . 

Jo 

where ymax is the outer boundary in the wall normal direction. The second was simply the wall pressure 
at a fixed streamwise location, Xi 

Ip = p(0,Xi)p*(0,Xi) . 

Note that both discrete cost functionals Ijj and Ip are essentially “delta functions”, being defined at only 
one streamwise gridpoint, which violates the necessary smoothness constraints needed during the derivation 
of the continuous adjoint. Qualitatively similar results were obtained when smooth cost functions were used 
to replace the the point functionals. Comparison of the discrete adjoint from grids with different resolutions 
requires a normalization step where the discrete adjoint is equated to its continuous counterpart by scaling 
it with a grid Jacobian. 


IV. Theoretical Receptivity 

The basic framework for predicting the roughness induced excitation of stationary crossflow modes is 
already available, at least for cases where the roughness amplitude is sufficiently small and the wavenum- 
ber spectrum of the surface height perturbations overlaps with the length scales of the crossflow vortex 
instabilities. Based on an extension of the Goldstein 6 -Ruban 7 theory for localized receptivity, the perturba- 
tion quantities associated with an arbitrary crossflow mode of spanwise wavenumber [3 can, in general, be 
expressed in the form: 10 

q(x,y,z-,/3) = C(x,p)4>(x,y)exp(i / ad£ + /3z) (12) 

Jx 0 
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where (#, y , z) denote the surface tangential, surface normal, and spanwise coordinates within an orthogonal, 
body fitted coordinate system; a(x; /3) is the ^-aligned wavenumber component of the crossflow mode of 
interest; and (j)(x, y ) represents an appropriately normalized local mode shape at station x. The cumulative 
coupling coefficient C(x ) then denotes the effective amplitude of the crossflow mode at x = xq] it is 
a function of both the intrinsic receptivity characteristics of the underlying boundary layer flow and the 
roughness geometry F(x,z)\ 


C(x) = f F(£, (3) A(£, /?) exp [-*0 P (£, /?)] d£ 

(13) 

J Xq 


i 

3 

ll 

3 

Tn 

3 

c^rv 

(14) 


where the lower limit of integration xq corresponds to the upstream end of the roughness distribution; F(£, /3) 
corresponds to the local amplitude of the relevant (spanwise) Fourier harmonic of the roughness height 
distribution F(x,z), while A characterizes the efficiency of the local receptivity to surface roughness, 
i.e. the local Green’s function with the response measured at the forcing location £. The combination of 
terms: A (£,/?) exp [— i0(£, /?)], is essentially an adjoint function that characterizes the impulse response, 
measured at the location x to forcing at location £. 

V. Results 

V.A. Incompressible Swept Hiemenz Flow 

The swept Hiemenz flow is an ideal configuration for concept study or code verification for crossflow instability 
in three-dimensional (3D) boundary layers. The basic flow has a linearly varying x profile, given by self- 
similar boundary-layer equations, and a uniform spanwise velocity component in the free- stream. The 
predictions of localized receptivity for this flow were presented in Refs . 15 and 19 using non-parallel receptivity 
theory plus linearized Navier- Stokes, and a linearized PSE framework, respectively. It was shown that neglect 
of mean-flow non-parallelism led to considerable over-prediction of the efficiency function magnitude for 
receptivity due to short scale variations in surface geometry or surface suction velocity, but the results based 
on linearized PSE were in excellent agreement with those based on linearized Navier- Stokes equations. See 
references 15,19,41,42 for a more detailed discussion of swept Hiemenz flow. 

The receptivity of the swept Hiemenz flow is considered first, beginning with a code-on-code (PSE/HLNS) 
validation exercise. The flow conditions are Re = 400 and u ; = 0 ; [3 = 0.4, based on the following 
nondimensionalizat ions : 

= S*X* ; S* = dU*/dX*\ x , = o ; l* = y/(u*/S) ; 

X = X*/r ; Re = W^l* /v* ; U x = U^/W^ = X/Re (15) 

(Here, Re denotes the Reynolds number based on the spanwise velocity component along the attachment 
line and the length scale based on the velocity gradient in the chordwise direction.) Figure 1(a) compares the 
growth rate a u of an x velocity disturbance, as computed using the LASTRAC 25 PSE and the present HLNS 
code. The growth rates are indistinguishable for this problem, verifying the parabolizing approximation in 
the PSE approach for problems of this type. Figure 1(b) compares the efficiency function associated with 
wall suction. The results agree well with those presented in references 15, 19. 

V.A.l. The Adjoint Solution 

Figure 2 (a) shows the Hiemenz flow 43 component of the adjoint solution for the HLNS equations. The 
cost functional is the magnitude of the velocity perturbation (Ijj) at x = 400, a location sufficiently 
downstream for the flow perturbations to be completely dominated by the least stable crossflow instability 
excited upstream. Plotted as a function of the x distance are the real and imaginary parts of the sensitivities 
as well as their magnitude envelope. Note that the data are plotted on a logarithmic- linear scale. The 
sensitivities of the u- and w- adjoint variables have a similar behavior. The most favorable location for a 
perturbation in the normal wall velocity is at approximately x — 160 — 170 [i.e. close to the neutral location: 
(x = 162) for crossflow instability with [3 = O.4.] The streamwise wavelength of at this location is 
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approximately 7.85. Also note that the wavelength of \I/ V increases linearly with x , as does the velocity [To© 
[see equation (15)]. 

Figure 2(b) shows the results of a grid refinement study for the Hiemenz flow, HLNS adjoint variable 
\I/ V . Note that the data are plotted on a logarithmic- logarithmic scale. Four grids are used in this study, 
with x resolution ranging from 1,920 — 5, 120 collocation points. All grids were distributed linearly in the 
x coordinate, to ensure a fixed number of points per wavelength. The normal resolution is held fixed at 55 
Chebyshev collocation points in all cases. The data collapses when the adjoint variables are scaled by the 
local grid metric (in this case, the inverse streamwise grid spacing), which allows each grid to be compressed 
independently in the vicinity of roughness elements. Note the rapid convergence of the adjoint variable in the 
upstream portion of the flow, the region of most relevance to receptivity prediction. The grid convergence 
is slowest for locations between the location used to define the cost functional (x = 400) and the beginning 
of the buffer domain {x = 450), i.e. the region of least relevance to receptivity prediction. The resolution 
of the adjoint variable is far less demanding in the vicinity of the localized roughness elements, x = 160, 
but requires greater resolution near the cost functional. In short, the grid resolution requirements needed to 
resolve the primal and adjoint variables can be very different. 

Figure 2(b) also shows the results for a cost functional placed at x — 410. The upstream results are 
indistinguishable (after appropriate scaling) from those obtained using x = 400, while the solutions near 
the functional are quite different. Thus, locating the cost functional at a position far downstream from the 
neutral point (in flow of this kind) is sufficient to obtain a meaningful adjoint solution. 

Surface roughness is an efficient mechanism for exciting crossflow disturbances if its spectral component 
at (or near) the disturbance wavelength is large. 11, 15,20,42 To this en q^ q f s desirable to study the impact of 
roughness elements of various sizes on a fixed wavelength disturbance, or conversely, to study the impact of 
a fixed size roughness element on disturbances of various wavelengths. The adjoint HLNS solution procedure 
provides an extremely efficient means to accomplish these parametric studies; the first step of which is 
obtaining a high-quality adjoint solution for all the disturbance wavelengths in question. Figure 3 shows the 
HLNS primal and adjoint solutions for the swept Hiemenz flow for the conditions Re = 400, a; = 0 for a 
decade variation in the spanwise parameter j3 = 0.1, 0.2, 0.4, 0.6, 0.8, 1.0. Figure 3(a) shows the influence of 
the parameter f3 on disturbance growth and/or decay. In all cases, a suction hole is located on the surface 
at the spatial location x = 200. Plotted on a linear- logarithmic scale is the variation in receptivity as a 
function of the forcing location, measured as a function of the magnitude of wall pressure fluctuation. The 
/ 3 = 0.4 disturbance has the largest amplitude at the outflow plane while that of f3 = 1.0 is the smallest. 
All disturbances decay super-exponentially ahead of the wall disturbance. Figure 3(b) shows the adjoint 
solution as a function of the parameter (3. Note that the adjoint solutions are qualitatively consistent 
with the growth rates shown in Fig. 3(a). 

V.A.2. Spanwise Periodic Roughness Arrays: DREs 

The definition of a spanwise periodic array of roughness elements has at least four independent parameters 
that can be optimized to accomplish a desired effect; 1) the roughness size and shape (e.g. cylindrical with 
radius r), 2) the roughness spacing (i.e. wavelength in spanwise direction), 3) the chordwise position on 
the airfoil, and 4) the height of the element. The first two are closely related (even dependent) and are 
used to change the relative mix of various harmonics in the boundary layer, while the third positions the 
disturbance at a desirable location (e.g. the neutral point in the stability curve). The fourth provides the 
amplitude of the disturbance. (Clearly the optimality of a particular DRE also depends on all the flow 
parameters Rc^cxaoA, • • • as well as the the spanwise mode (/?) in question.) We begin by optimizing the 
diameter/shape, spacing and position of a spanwise DRE array used to excite swept Hiemenz flow. 

The shape of each element is assumed to be a circular cylinder, with the corners of the cylinder assumed 
to be sharp; rounding them had little impact. The radius of the circle was used as the first parameter in the 
study. The spanwise spacing for this exercise was chosen to coincide with the wavelength of the spanwise 
mode. Next, given a [3 (and A = 27r//3), an FFT was used at each x position to expand the spanwise 
shape into its Fourier components, whereby the amplitude at the wavelength A was identified. Note that this 
amplitude varies with x position, depending on the relative position on the element. For example, assume 
the radius to be half the spanwise period: 2 r = A. The spanwise shape in the middle of this element xq is 
a constant because the elements exactly touch in this case. All off-center spanwise shapes are square waves 
of shrinking breadth (with period A), until the positions x = xo db r are reached, and are zero thereafter. 

Figure 4 shows the effect of DRE radius and location on the amplitude of the crossflow mode corresponding 
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to the fundamental harmonic of the array spacing. The flow conditions are Re = 400 and uj = 0. Results 
for three spanwise wavenumbers {(3 = 0.2, 0.4 0.8) are shown in Figs. 4(a), 4(b), and 4(c), respectively. 
All results were obtained by convolving the adjoint solution for each respective [3 with the DRE surface 
boundary condition. The adjoint cost functional used for this study was the surface pressure /p, and was 
located sufficiently downstream to not influence the results. The abscissa in each plot is the element radius 
normalized by the spanwise wavelength A = / (3 being considered. The ordinate is the amplitude of the 

surface pressure at the objective location. (Note the scale changes for each (3). Five x positions were tested 
for each spanwise (3 , centered approximately at the neutral point. Observations include 

• The peak sensitivity coincides in all cases at the neutral point, with a modest decrease in sensitivity 
for upstream or downstream locations. This result follows immediately from the breadth of the adjoint 
solution [see Fig. 3(b)] in the vicinity of the neutral point. 

• Cyclic peaks and attenuation’s occur as the radius of the DRE is varied from 0 < r < A. 

V.A.3. Deterministic Nonlocal Roughness: Wavy Wall 

As a building block for distributed receptivity predictions, we extend the results in Ref. 19 to surface waviness 
of the form F(x,z) = exp [i(a w x + f3z)\ using both the linearized PSE and HLNS frameworks. (Only the 
HLNS results are presented herein, as the results obtained by both approaches are nearly indistinguishable; 
a further verification of the adjoint PSE predictions.) As discussed in Ref., 43 predictions for such wavy wall 
configurations can be easily applied to estimate the receptivity due to a spatially homogeneous roughness 
distribution. 

The x evolution of the effective coupling coefficients for 13 values of a w is shown in Fig. 5, calculated by 
convolving the Hiemenz wall adjoint solution with the individual surface waviness boundary condition. Each 
profile a wj , j = 1, 13 is plotted using the detuning parameter d = (a w - a neutra i) / a neutra i where a neutra i 
is the local wavenumber at the lower branch neutral point. This definition of d normalizes the individual 
difference in wavenumber a w — ot ne utrai against a neutra i. Figures 5(a) and 5(b) display the spatial evolution 
of the effective coupling coefficient of the positively and negatively detuned disturbances, respectively, for 
the spanwise wavenumber f3 = 0.4. Figures 5(c) displays the coupling coefficient at the spatial location 
x = 400 as a function of the detuning parameter d, for four spanwise wavenumbers: f3 = 0.2, 0.4, 0.6, 0.8. 

Figures 5(a) and 5(b) show that the receptivity to wavy walls is effectively confined to a relatively 
narrow region, centered on the location where the wavenumber a w of the surface height distribution most 
closely matches the real part of the local instability wavenumber in the x direction. Because chordwise 
wavenumber of the crossflow mode decreases with x, the wavenumber synchronization at smaller a w occurs 
farther downstream, consistent with the progressively downstream movement of the region of rapid rise in 
C as the detuning parameter d becomes more negative [see Fig. 5(b)]. The variation in the overall coupling 
coefficient (i.e., the value of the downstream asymptote of C) with respect to d is roughly symmetric about 
d — 0, although the actual peak is slightly shifted to d = —0.10. This symmetric dependence on the detuning 
d is similar to the case of acoustic excitation of Tollmien-Schlichting (TS) waves due to small amplitude 
waviness imposed on a flat surface. 43 The striking difference between the two cases, however, corresponds 
to the range of the detuning parameter over which relatively significant receptivity occurs. Whereas the 
coupling coefficient becomes smaller than 10 percent of the peak value for |d| > 0.10 for TS excitation over 
a flat plate, the swept Hiemenz flow appears to be receptive to surface waviness for detuning parameters in 
excess of \d\ > 0.5. 

Figure 5(c) shows the overall effect of detuning at several spanwise wavenumbers (3. Here, the downstream 
asymptote of the effective coupling coefficient is plotted as a function of the detuning parameter d. Unlike 
the (3 = 0.4 case, detuning lacks symmetry for other values of (3. Indeed, the decay in C for positive detuning 
is significantly slower than for negative detuning. 

Table 1 compares the receptivity of local DREs and distributed wavy wall surface boundary conditions, 
for three values of the parameter f3. The actual numbers in the table correspond to 1) the maximum response 
to local DREs of any size or surface location [i.e. the maximum amplitude derived from Figs. 4(a), 4(b), 
and 4(c)] and 2) the maximum response to a distributed wavy wall for any value of detuning parameter d 
[i.e. the maximum amplitude derived from Fig. 5(c)]. Recall that the adjoint cost functional used for both 
studies is the surface pressure Ip, which implies that this comparison relates the surface pressure at x = 400: 
a location far downstream from the region of peak sensitivity. Clearly, the distributed wavy wall is much 
more effective at exciting the growth of disturbances; over an order of magnitude in all cases. 
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P = 0.2 
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II 

o 

4^ 

P = 0.6 

DREs 

6e-6 

le-4 

7e-6 

Wavy Wall 

le-4 

4e-3 

3e-4 


Table 1. Comparison of the amplitude of the pressure disturbance at x = 400 achieved using DREs or wavy wall surface 
profiles . 


V.B. SWIFT Configuration 

A long term goal is to predict the initial crossflow amplitudes via a reduced order specification of the roughness 
geometry, using only a small number of amplitude and texture parameters for the surface roughness. It is 
envisioned that the study will help to identify the most relevant roughness specification parameters and hence, 
contribute to the potential development of an appropriate database pertaining to the natural roughness of 
typical aerodynamic surfaces. 17 Here, we report the progress that has been made to date. 

The study begins with a general description of the stability characteristics of the SWIFT configuration. 
Next, the adjoint solution is used to identify the highly sensitive regions near the nose, followed by a 
parametric study that identifies the optimal radius and placement of DREs. Finally, experimentally obtained 
surface roughness profiles are used to assess the receptivity characteristics. Receptivity predictions for the 
suction surface of the SWIFT configuration are presented. The flow conditions correspond to Test Point 
27, with —4.69° angle of attack and Re c = 7.155 x 10 6 . The selected spanwise station corresponds to an 
inboard section at 13 inches from the root (i.e., mounting) location of the test article. 

Figure 6 shows the evolution of the crossflow disturbance amplitude in the chordwise direction x. The 
abscissa is the nondimensional distance from the stagnation streamline, while the ordinate is the amplitude 
of the energy 

nymax 

E{x) = 1/2 / [u(y,x)u*(y,x) + v(y,x)v*(y,x) + w(y,x)w*(y,x)]dy 

Jo 

associated with a crossflow mode with a spanwise wavelength of 4.5mm. Figure 6(a) shows the disturbance 
growth on the scale of the chord, while Fig. 6(b) shows a magnified view of the nose region. Each plot 
compares the solutions obtained with the HLNS solver on grids of different density, and that obtained from 
the linear PSE solver LASTRAC. 19 Excellent agreement is achieved throughout most of the domain, until 
the buffer domain region (0.55 < x ) is encountered in the HLNS formulation. The PSE calculation begins 
near x ~ 15mm, close to the neutral point of the crossflow mode. The HLNS solution corresponds to the 
forced response of the boundary layer to a spanwise periodic array of DREs near x ~ 10mm. Because of the 
linearity of both solutions, the PSE solution is scaled to match the HLNS solution at x = 300mm. 

V.B.l. The Adjoint Solution 

Figure 7(a)-(b) shows the solutions to the adjoint HLNS equations for the surface variables T u , 4/ v , The 
objective function used for these solutions was the surface pressure Ip at the point x = 0.5 chord. The real 
and imaginary components and the magnitude of each surface variable, e.g. [Re(Sk u )], [Jm(T n )], and |T U |, 
are plotted. Figure 7(b) is a coarser view of Fig. 7(a) used to compare the sensitivities and to forcing 
via the u and w wall components, respectively. The most sensitive portion on the surface of the airfoil, is 
located upstream of the x = 2% chord (i.e. x < 27mm) and then decreases rapidly with distance downstream 
of this peak location. The significance of receptivity within the leading edge region is consistent with both 
previous experiments 22 and theoretical predictions. 12, 16,17,19 Also note that the chordwise wavelength of 
the adjoint components decreases rapidly in this region. 

Figure 8 extends the results shown in Fig. 7 to include the influence of different spanwise wavelengths. 
Shown are the amplitudes of the adjoint variable over the first 100mm of the chordwise distance. Spanwise 
wavelengths in the range 2.25mm < A < 21.5mm are included in the study. Figure 8(a) shows the amplitude 
of the adjoint variable \^f v \ for all wavelengths, as well as Re[^/ V ] for the A = 2.25mm and A = 21.5mm 
cases. The wave nature of the chordwise wall sensitivity is clearly evident in the Re[^t v ] profiles shown for 
A = 2.25mm and A = 21.5mm. Note that the chordwise wavelength of the sensitivity, scales with the 
spanwise wavelength A near the nose and expands rapidly beyond approximately x = 10mm. 
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V.B.2. Spanwise Periodic Roughness Arrays: DREs 

Figure 9 shows convolutions of the adjoint variables shown in Fig. 7 against surface boundary conditions 
approximating DREs. The parameter investigation is equivalent to that performed for the swept Hiemenz 
flow (see Fig. 4); specifically, the shape of each DRE is assumed to be a circular disk of variable radius r, with 
spanwise spacing A = 2 tt//3. The position of the DRE array is varied over the interval 0 < x < 8 A, which 
spans the neutral point on the surface of the airfoil and most of the highly receptive region near the nose. 
The abscissa (linear) of Fig. 9 shows the radius of the DRE normalized by the spanwise wavelength A, while 
the ordinate (linear) shows the amplitude of the wall pressure at the mid-chord location. (All calculations 
were performed using the adjoint HLNS formulation, with the Ip cost functional situated mid-chord.) 

The spanwise oriented DRE array is highly effective at exciting the crossflow disturbance within the 
range x ~ (3 — 4)A which is 13.5 — 18mm from the attachment line. Note that this location is close to 
the neutral point identified by the PSE computation shown in Fig. 6 (b). The amplitude in this region [e.g. 
x = (2.0 — 4.0) A] is relatively uniform, but decreases rapidly outside the interval. The receptivity falls to 
near zero as the radius of the DRE approaches 2r = A. 

In general, the DRE results for the SWIFT airfoil are similar to those obtained for DREs subjected to 
the swept Hiemenz mean flow. Specifically, the peak amplitude coincides with the peak of the adjoint, and 
there is slow decay in peak amplitude for upstream or downstream locations. Subtractive cancellation occurs 
as the radius of the DRE traverses through resonant wavelengths. The region 13.5 — 40mm seems a likely 
candidate for a wavy wall surface study, and is the object of ongoing work. 

V.B.3. Random Surface Roughness: Experimental Data 

The final exercise is to establish the disturbance behavior resulting from physically realistic surfaces. Several 
different surface textures/finishes were applied to the SWIFT airfoil to study the impact of these textures 
on transition. The data were provided by Saric , 44 taken from an ongoing Texas A&M University rough- 
ness/transition experiment. The surface textures included regions with the following distinguishing features: 
1) a taped surface, 2) a polished surface, 3) a painted surface, and 4) a painted surface that was subsequently 
sanded. The physical characteristics of each texture were obtained experimentally using a stylus that sweeps 
linearly along the spanwise direction. The sweep length is 1.6 inches over which 5,000 datapoints were 
collected. In total, 13 distinct surface samples were collected from four texture groups. 

Table 2 presents the data used to quantify the 13 surface samples, and as expected each texture has its 
own characteristics. Cases 1—4, 5 — 9, 10 — 11 and 12 — 13 were obtained from the 1 ) taped, 2 ) polished, 3) 
painted and 4) painted/sanded categories, respectively. The metrics used to characterize the surface textures 
are described in the caption of the figure, e.g. = raw mean surface absolute deviations. Several 

observations about the surface data are worthy of note. 

• The painted surface has by far the roughest texture, while the polished surface is the smoothest. They 
differ in RMS surface deviation by over an order of magnitude. 

• There is considerable deviation within the raw data for each individual texture group. For example, 
the raw data reported in case 07 is significantly different from the other members of the polished group. 

• A spectral decomposition is needed to quantify the amplitude of each wavenumber component of the 
surface height distribution. Because the data are neither periodic nor smooth, this decomposition is 
problematic if Fourier techniques are used. Thus, a detrending step (described shortly) or equivalently 
a high pass filter, is necessary to correct for the large scale variations comparable to the length of each 
scan. 

• Detrending the data profoundly changes the values of some of the metrics used to characterize the 
surface texture. Apostori testing must be used to establish the efficacy of the detrending step. 

• Changes in raw maximum distance peak to peak: £3 and raw five point average maximum distance 
peak to peak: (4, as a result of detrending were not checked. 

Figure 10(a) shows an example of how the surface data are detrended. The data being detrended is 
reported in Table 2 as Case 3 , and belongs to the taped group. In Fig. 10(a), the raw data (shown in 
red) has significant long wavelength variation. Superimposed along this long wavelength profile are high 
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Raw Data 



Detrended 


Case 

|| hr aw Hi 

|| hr aw || 2 

(3 

(4 

\\hdet Hi 

\\hdet || 2 


01 

0.84 

0.99 

3.97 

3.66 

0.06 

0.08 

Taped 

02 

1.32 

1.55 

5.49 

5.36 

0.09 

0.11 

03 

1.29 

1.52 

5.44 

5.32 

0.08 

0.11 


04 

0.82 

1.00 

4.69 

3.85 

0.09 

0.12 


05 

0.12 

0.17 

1.53 

0.92 

0.08 

0.11 


06 

0.12 

0.15 

1.14 

0.86 

0.06 

0.08 

Polished 

07 

0.31 

0.36 

1.89 

1.18 

0.09 

0.11 


08 

0.15 

0.17 

1.10 

0.84 

0.07 

0.10 


09 

0.15 

0.18 

1.19 

0.74 

0.07 

0.08 

Painted 

10 

6.27 

8.06 

48.12 

42.22 

6.06 

7.86 

11 

4.37 

6.00 

49.16 

30.16 

4.35 

5.93 

Painted/ Sanded 

12 

1.08 

1.18 

16.25 

13.49 

0.58 

1.08 

13 

0.95 

1.78 

6.14 

3.29 

0.50 

0.95 


Table 2 . SWIFT experimental data: surface roughness statistics. The column designations are as follows; ||/i r a«;|li = 
Raw mean of surface absolute deviations; \\h raw ||'g = Raw root mean square of surface deviations; £3 = Raw maximum 
distance peak to peak; £4 = Raw five point average maximum distance peak to peak; ||/idet||i — Detrended mean of 
surface absolute deviations; || || 2 = Detrended root mean square of surface deviations. 


frequency surface oscillations which appear to be random. Detrending the data begins by first projecting the 
raw data onto a low-order polynomial using an L 2 minimization of the deviation. This allows the coefficients 
of a long wavelength polynomial (shown in green), to be determined. Subtracting the low-order polynomial 
from the original data produces a detrended dataset (shown in blue), that is approximately periodic. The 
polynomial order used to detrend all data in this study was 4; polynomials of degree higher than 4 decrease 
the energy of the modes of most interest, while those of lower degree were inadequate for the purpose of 
detrending. 

Figure 10(b) compares the Fourier spectrum of the original raw data (solid lines), with that of the 
detrended data (dashed lines). The four cases shown all originate from the polished surface group. Note 
that detrending the original data is highly effective at reducing the first three Fourier modes (A > 15mm), 
while leaving the high frequency character of all other modes. Furthermore, note that the Fourier modes 
corresponding to the most unstable mode of the crossflow instability (4.5mm) remain virtually unaffected 
by the detrending step. 

Figures 11 and 12 show the original raw and detrended data for the polished cases and the painted and 
painted/sanded cases. Note that detrending the painted data has little effect on the apparent shape of the 
surface profiles. This is consistent with the data in Table 2 which indicates that the raw surface measures: 
|| Kaw 111, \\Kaw\\ 2 i closely match the detrended measures: \\h det ^ lf ||hdet|| 2 - 

V.B.4- Random Surface Roughness: Disturbance Growth 

To determine the impact of the different textures on crossflow disturbance amplitude, the near leading edge 
surface of the SWIFT airfoil is fitted with each of the 13 experimentally measured roughness patches. The 
region 0 < x < 0.02 % chord (0 < x < 27 mm) was shown previously (see Figs. 7 and 9 herein), to be the 
most efficient at exciting the crossflow disturbances. The patches extend approximately 40mm chordwise, 
allowing each to completely cover the sensitive nose region. 

An assumption is necessary in the chordwise direction to fully specify the surface profile. For brevity, 
roughness characteristics along the chordwise direction are assumed to be the same as those in the spanwise 
direction, and that the airfoil surface is composed of a tensor-product roughness patch for each of the 13 
surface roughness patches. 

Using the surface roughness patches as boundary conditions, the adjoint HLNS methodology is then used 
to predict the crossflow disturbance amplitude for each texture. Convolving each surface boundary condition 
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with the adjoint solution, \E^, yields the crossflow disturbance amplitude at the location of the specified cost 
functional. The Ijj cost functional is used exclusively in this study and is always located at the mid-chord 
position. 

Figures 13 and 14 present an extensive study of disturbance sensitivity to surface roughness textures. 
The abscissa (linear) of each figure covers the spanwise wavenumber interval 100 < (3 < 5000, (i.e. 
86mm > A* > 1.7mm; lengths nondimensionalized by 4.5’ chord) approximately centered about the highly 
unstable [3 = 1915 wavenumber (i.e. the 4.5mm spanwise wavelength). The ordinate (logarithmic) measures 
the disturbance amplitude for each test case. 

Each of the two figures contains 611 symbols; 13 symbols (each corresponding to an individual surface 
texture patch) at 47 distinct wavenumbers /3. The symbols for the 13 surface texture patches form a vertical 
distribution at each wavenumber /3. Recall that the 13 surface textures were divided into four groups: taped, 
polished, painted, and painted/sanded. Distinct symbols (and colors) are used for each surface texture group. 
For example, the green inverted triangles denote the disturbance amplitude resulting from the two painted 
surface patches at each of the 47 wavenumbers. 

Each of the 611 symbols from Figs. 13 and 14 is accompanied by an overbar - . The over-bars quantify 
the 1) ensemble mean disturbance amplitude (with 128 realizations) and 2) RMS deviation in response to 
one surface roughness patch. (The bottom half of the RMS deviation bar is omitted for reasons of clarity.) 
The ensemble of roughness distributions is constructed as follows. First, note that the phase of the data in 
each roughness patch is completely arbitrary because of periodicity (an assumption shown to be essentially 
valid in the detrended case). Thus, an arbitrary number of the disturbance responses can be obtained from 
the data in a single roughness patch. For example, assume a roughness patch extends over the surface 
interval 0 < x < 40 mm, but that the data are translated periodically N times, such that the virtual origin 
xo (the first data point) satisfies the relation xq = (k — 1) 40rara/7V, k = 1,7V. Indeed, the shift could even 
be random. Using this periodic translation of boundary data, an ensemble mean and standard deviation 
in disturbance amplitude can be determined for each surface patch. Numerical experimentation determined 
that the mean and RMS are converged to three significant digits using an ensemble of 128 samples. 0 

Figures 13 and 14 are obtained using the raw and detrended data, respectively. Note that there is little 
difference in disturbance amplitudes for wavenumbers in the interval (3 > 1300(^ 6.5mm) and no difference 
for [3 > 1915(« 4.5mm). The differences in receptivity predictions without and with detrending of surface 
height distributions, are consistent with the associated variation in the Fourier spectra of the surface height 
distribution. Specifically, the effect of detrending is to decrease the amplitudes of the long wavelength modes 
(e.g. A > 15mm), which are outside the wavelength band that is relevant for transition as well as transition 
control via DREs. 

The following observations summarize this study on the effects of random surface roughness on the 
strength of crossflow disturbance excitation. 

• The most striking result is the wide range of crossflow amplitudes resulting from the four surface 
textures. Approximately two orders of magnitude variation in disturbance amplitude is observed when 
comparing the painted and polished surfaces. The differences in predicted disturbance amplitudes for 
the painted/sanded leading edge and the polished edge are smaller, but still approximately one order 
of magnitude. 

• In addition to the surface metrics, \\h raw \\ k ,k = 1,2 ; £/c, k = 3,4 and \\hd e t\\k^ = 1? 2, an accurate 
Fourier decomposition of the surface roughness is needed to predict the disturbance response at all 
wavenumbers. Non-periodicity and aliasing can lead to erroneous component amplitudes. 

• The range of mean disturbance amplitudes (max - min) arising from the individual members of a texture 
group, greatly exceeds the RMS variation in the individual ensemble averages within that group. This 
effect could be greatly mitigated by surface samples of significantly greater length. Furthermore, the 
in-group deviation depends strongly on the disturbance wavenumber. 

The effect of surface quality on transition has been documented in Ref. 1, wherein transition was noted 
to move from somewhere between 25 to 30 percent chord to 80 percent chord when the painted leading 
edge was polished to reduce the roughness amplitude from 1.0 fi m RMS (3.8 /im average peak-to-peak) to 
0.3 fi m RMS (2.2 /i m average peak-to-peak). Had linear stability based N-factor methods been applied to 

c Figures 13 and 14 are each the result of 128 X 13 X 47 = 78, 208 individual convolutions of surface data with the 47 adjoint 
solutions. The adjoint HLNS capability is enabling for problems of this complexity. 
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predict the onset of transition in this case, the correlating N-factor would have increased from 8 to 16. The 
substantially smaller initial crossflow amplitudes on the polished leading edge (Fig. 14) are qualitatively 
consistent with the large delay in measured transition location. 

Transition prediction via N-factor correlations for high-frequency secondary instabilities of the stationary 
crossflow modes was considered in Ref. 24 in the context of modeling transition delay via discrete rough- 
ness elements (DREs). By applying the same methodology in the absence of any DREs, transition onset 
location was predicted to shift downstream by approximately 15 percent of the chord length for each order- 
of-magnitude reduction in the initial amplitude of the dominant 4.5mm stationary crossflow mode. The 
predicted shift in transition location could be considerably larger if only one of the two dominant types of 
secondary instability modes (namely, the family of secondary instability modes that is associated with span- 
wise gradients of the basic state, as against wall-normal gradients) were to be realized in practice. Because 
the predicted impact of the reduced crossflow receptivity on the polished leading edge appears to be smaller 
than that measured during the SWIFT experiments, further work is necessary to fully explain the observed 
transition behavior. 


VI. Summary 

This paper extends the prediction capability for roughness induced crossflow receptivity to realistic surface 
roughness distributions, to pave the way for integrated, high fidelity analysis methods for crossflow transition, 
in general, and DRE based laminar flow control, in particular. An illustrative set of results relevant to a 
unique experiment involving a swept wing configuration were presented above. 

The monumental scope of the problems at hand (e.g. quantification of flow response to random surface 
roughness via Monte-Carlo simulation of canonical and realistic geometries) made it expedient to develop an 
adjoint HLNS capability, to complement an existing HLNS solver. A discrete adjoint approach was adopted 
because of complexities resulting from the artificial boundary conditions used at the outflow plane of the 
HLNS solve and the fact that the discrete adjoint is an exact adjoint to the forward problem in question. 
The new adjoint formulation was exercised on a set of well known crossflow stability problems including the 
swept Hiemenz flow. 

The new adjoint formulation was used to tackle the problem of distributed receptivity: excitation of 
stationary crossflow modes via a wavy wall F(x,z) = exp [i(a w x + f3z)] was investigated in the context 
of swept Hiemenz flow. The receptivity to wavy walls is effectively confined to a relatively narrow region 
centered on the location where the wavenumber a w of the surface height distribution most closely matches 
the real part of the local instability wavenumber along the chordwise x direction. The variation in the 
overall coupling coefficient with respect to a detuning parameter d = (a w — ctneutrai) / ^neutral, is roughly 
symmetric about d — 0 (although the actual peak is slightly shifted to d = —0.10), similar to the case of 
acoustic excitation of TS waves due to waviness over a flat surface. 43 However, the striking difference between 
the two cases corresponds to the range of detuning parameter over which relatively significant receptivity 
occurs. Whereas the coupling coefficient becomes smaller than 10 percent of the peak value for \d\ > 0.10 
for TS excitation over a flat plate, the swept Hiemenz flow appears to be receptive to surface waviness for 
detuning parameters in excess of \d\ > 0.5. 

Finally, the adjoint formulation was used to study the SWIFT experiment. One long term goal of this 
study is to enable an estimation of initial crossflow amplitudes via a reduced order specification of the rough- 
ness geometry, in the form of a small number of amplitude and texture parameters for the surface roughness. 
The study began with a general description of the stability characteristics of the SWIFT configuration, and 
confirmed that the most sensitive region is near the neutral point. DREs located in this region have the 
most impact. Next, the adjoint solution was used study the impact of different surface roughness patches on 
disturbance growth of unstable crossflow modes. An exhaustive exercise for 47 distinct spanwise wavelengths 
resulted in the following observations: 1) two orders of magnitude variation in disturbance amplitude are 
observed when comparing the different surface textures (painted vs. polished surfaces). This is qualitatively 
consistent with the significant shift in predicted transition location when the initial crossflow amplitude is 
reduced by one to two orders of magnitude. 2) Accurate Fourier decomposition of the surface roughness 
is needed to predict the disturbance response at all wavenumbers. 3) Significant variation in disturbance 
amplitudes appears within each of the four texture groups, and is greater than the RMS deviation predicted 
within each of the 13 individual surface samples. Much work still remains to establish a reliable reduced set 
of surface metrics. 
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Much work is still left to be done regarding the quantification of surface roughness. The results presented 
herein are encouraging in suggesting that a low dimensional surface parametrization could eventually be 
identified. Ultimately, modeling random, spatially distributed small magnitude roughness on this and other 
canonical configurations will aid in the development of improved metrics (or improved interpretation of 
existing set of metrics) for the specification of natural surface roughness (e.g., RMS, peak-to-peak, spectral 
probability density function model); the end result being a more reliable set of assessment /design tools for 
DRE configurations and fast estimates of initial crossflow amplitudes for nonlinear PSE calculations with or 
without DREs. 
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Figure 1. Comparison of the PSE and HLNS formulations on swept Hiemenz flow, Re = 400, (3 = 0.4 with suction 
surface boundary conditions. 
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(a) Wall sensitivity to normal perturbations. (b) Grid convergence of discrete adjoint. 

Figure 2. Adjoint solution for the swept Hiemenz flow for conditions Re = 400, (3 = 0.4. Figure (a) shows the surface 
sensitivities of the real, imaginary and magnitude of 4> v ; the cost functional Ijj is situated at x = 400. Figure (b) shows 
a grid refinement study of the adjoint surface variable 4> v , as well as the influence of the location of Ijj on 4^. 




(a) Primal HLNS: Wall pressure for various spanwise / 3 . (b) Adjoint HLNS: 4 ^ v - sensitivity for various / 3 . 

Figure 3. HLNS primal and adjoint solutions for the swept Hiemenz flow: (Re = 400, uj = 0), comparing the influence 
of (3 on disturbance growth. Figure (a) shows the (3 sensitivity on growth and/or decay of disturbances initiated by wall 
suction located at x = 200. Note that the buffer domain begins at x = 950. Figure (b) shows the adjoint solution 4' ?; as 
a function of the same (3 parameters. 
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(a) HLNS: (3 = 0.2. 


(b) HLNS: (3 = 0.4. 


(c) HLNS: (3 = 0.8. 


Figure 4. HLNS: Radius study performed on the swept Hiemenz flow for conditions Re = 400, u> = 0, showing sensitivity 
of radius on spanwise wavelength A. 



(a) HLNS: Positive Detuning. 



(b) HLNS: Negative Detuning. 



Figure 5. HLNS: Detuning study performed on the swept Hiemenz flow for conditions Re = 400, (3 = 0.4 are shown in 
Figs, (a) and (b). The detuning parameter is defined as d = (ot w —a neu tral)/ a neutral where cx w is the chordwise wavelength 
of the surface disturbance. Figure (c) shows a detuning study for the spanwise wavenumbers (3 = 0.2, 0.4, 0.6, 0.8. 




X (Nondim by Chord) 


(a) Disturbance Energy Growth: Full 


(b) Disturbance Energy Growth: Zoom 


Figure 6. Growth of the disturbance energy with respect to chordwise position x , comparing linear PSE and HLNS. 
Flow conditions are Re = 7.15 X 10 6 , (3 = 1915 (4.5mm wavelength) and u> = 0. Two grid resolutions are used for the HLNS 
simulations to assess convergence. 


18 of 23 


American Institute of Aeronautics and Astronautics 


— -A— - Re(u) 
A lm(u) 



10 ° 10 1 10 2 

X (mm) 


(a) Adjoint HLNS: ^ u solution (b) Adjoint HLNS: ^ v solution 



250000 


250000 


10 ° 10 1 10 2 

X (mm) 


— -A— - Re(u) 

A lm(u) 

Re(v) 

0 lm(v) 

A -|u| 

H -|v| 

— e— -iwi 

A +|u| 

■B +M 

-0 +|W| 


Figure 7. Adjoint HLNS solutions for the SWIFT airfoil flow for conditions Re = 7.15 X 10 6 , (3 = 0.0314 (4.5mm 
wavelength). The region of maximum sensitivity is approximately 10mm and decreases rapidly downstream. The ^f v is 
the most sensitive component. 
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Figure 8. The |4A| solution is presented for 7 spanwise wavelengths, centered about 4.5mm. The real component I?e[4A] 
is also included for the shortest and longest wavelengths, to show the dependency of (chordwise) wavelength on position 
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Figure 9. The sensitivity to DRE radius and location is examined for the SWdFT airfoil flow for conditions Re = 7.15xl0 6 , 
/ 3 = 1915 (4.5 mm wavelength). The chordwise region of maximum sensitivity is about 3 wavelengths, and coincides 
approximately with the neutral location. 




Figure 10. Figure (a) demonstrates the process of detrending on a taped a dataset. First the raw data are fitted 
with a low-order polynomial, which is then subtracted from the original data to yield a detrended dataset. Figure (b) 
shows the impact of detrending on the Fourier coefficients of a polished dataset. Detrending primarily affects the long 
wavelength modes. 
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(a) Raw Polished Surface Data 


(b) Detrended Polished Surface Data 


Figure 11. Line data for the five surfaces in the polished group. Figure (a) shows the raw surface data, while Fig. (b) 
shows the data after detrending. 



(a) Raw Painted Surface Data 



(b) Detrended Painted Surface Data 


Figure 12. Line data for the two surfaces in the painted group, and three surfaces of the painted/sanded group. Figure 
(a) shows the raw surface data, while Fig. (b) shows the data after detrending. 
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WaveNumber (P) 

Figure 13. The Adjoint HLNS formulation is used to simulate the crossflow disturbance response to four classes of 
surface roughness: Taped, polished, painted and painted/sanded. Re = 7.15 X 10 6 , uj = 0. The adjoint cost functional is 
placed at 50% chord. Raw surface data are used in the study. 
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Figure 14. The Adjoint HLNS formulation is used to simulate the crossflow disturbance response to four classes of 
surface roughness: Taped, polished, painted and painted/sanded. Re = 7.15 X 10 6 , uj = 0. The adjoint cost functional is 
placed at 50% chord. Detrended surface data are used in the study. 
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